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Abstract 

We present 0{al) corrections to deep inelastic scattering amplitudes on massive quarks 
obtained within the scheme of Aivazis, Collins, Olness and Tung (ACOT). After iden- 
tifying the correct subtraction term the convergence of these contributions towards the 
analogous coefficient functions for massless quarks, obtained within the modified minimal 
subtraction scheme (MS), is demonstrated. Furthermore, the quantitative relevance of 
the contributions to neutral current (NC) and charged current (CC) structure functions 
is investigated for several choices of the factorization scale. 



1 Introduction 



Leptoproduction of heavy quarks has become a subject of major interest in QCD phe- 
nomenology both for experimental and theoretical reasons. Heavy quark contributions 
are an important component of measured neutral current (NC) [|l], and charged current 
(CC) deep inelastic (DI) structure functions at lower values of Bjorken-x, accessible 
to present experiments. Charm tagging in NC and CC deep inelastic scattering (DIS) 
offers the possibility to pin down the nucleon's gluon and strange sea p, ||, density, 
respectively, both of which are nearly unconstrained by global fits to inclusive DI data. 
Theoretically it is challenging to understand the production mechanism of heavy quarks 
within perturbative QCD. The cleanest and most predictive method |^ of calculating 
heavy quark contributions to structure functions seems to be fixed order perturbation 
theory (FOPT) where heavy quarks are produced exclusively by operators built from 
light quarks (u,d,s) and gluons (g) and no initial state heavy quark lines show up in 
any Feynman diagram. Heavy quarks produced via FOPT are therefore also called 'ex- 
trinsic' since no contractions of heavy quark operators with the nucleon wavefunction 
are considered (which in turn would be characteristic for 'intrinsic' heavy quarks). Be- 
sides FOPT much effort has been put into formulating variable flavor number schemes 
(VFNS's) |11|, |I2], which aim at resumming the quasi-collinear logs [ln((5^/m^); 

Q and m being the virtuality of the mediated gauge boson and the heavy quark mass, 
respectively] arising at any order in FOPT. All these schemes have in common that ex- 
trinsic FOPT induces the boundary condition |T^, |TB|] g(x, = rn?) = + 0{a1) for an 
intrinsic heavy quark density, which then undergoes massless renormalization group (RG) 
evolution. Apart from their theoretical formulation VFNS's have to be well understood 
phenomenologically for a comparison with FOPT and with heavy quark tagged DI data. 
We will concentrate here on the scheme developed by Aivazis, Collins, Olness and Tung 
(ACOT) 0, 0. In the ACOT scheme full dependence on the heavy quark mass is kept in 
graphs containing heavy quark lines. This gives rise to the above mentioned quasi-collinear 
logs as well as to power suppressed terms of 0[{rn? /Q'^Y]. While the latter give mass 
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corrections to the massless, dimensionally regularized, standard coefficient functions (e.g. 
in the MS scheme), the former are removed by numerical subtraction since the coUinear 
region of phase space is already contained in the RG evolution of the heavy quark density. 
Up to now explicit expressions in this scheme exist for DIS on a heavy quark at 0{a^) 
as well as for the production of heavy quarks via virtual boson gluon fusion (GF) at 
0{al) [|^. In section 2 we will give expressions which complete the scheme up to 0{al) 



and calculate DIS on a heavy quark at first order in the strong coupling, i.e. B*Qi — > 
(inch virtual corrections to B*Qi Q2) with general couplings of the virtual boson B* 
to the heavy quarks, keeping all dependence on the masses mi 2 of the quarks Qi,2- It 
is unclear whether (heavy) quark scattering (QS) and GF at 0{al) should be consid- 
ered on the same level in the perturbation series. Due to its extrinsic prehistory QS*^^^ 
(bracketed upper indices count power^ of a^) includes a collinear subgraph of GF'-^\ e.g. 
7*c —>■ eg contains the part of '~f*g —>■ ccg, where the gluon splits into an almost on-shell 
cc pair. Therefore QS at 0{al) can be considered on the level of GF at 0{al). On the 
other hand the standard counting for light quarks is in powers of ag and heavy quarks 
should fit in. We therefore suggest that the contributions obtained in section 2 should be 
included in complete experimental and theoretical NLO-analyses which make use of the 
ACOT scheme. Theoretically the inclusion is required for a complete renormalization of 
the heavy quark density at 0{al). However, we leave an ultimate decision on that point 
to numerical relevance and present numerical results in section 3. Not surprisingly they 
will depend crucially on the exact process considered (e.g. NC or CC) and the choice of 
the factorization scale. Finally, in section 4 we draw our conclusions. Appendices A and 
B outline the calculation of real gluon emission and virtual corrections, respectively, and 
some longish formulae are presented in Appendix C. 



^ For the reasons given here we refrain in most cases from using the standard terminology of 'leading' 
and 'next-to- leading' contributions and count explicit powers of as- 
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2 Heavy quark contributions to structure functions 



In this section we will present all contributions to heavy quark structure functions up 
to 0{al). They are presented analytically in their fully massive form together with the 
relevant numerical subtraction terms which are needed to remove the collinear divergences 
in the high limit. Section 2.1 and 2.3 contain no new results and are only included 
for completeness. In section 2.2 we present our results for the massive analogue of the 
massless MS coefficient functions Cf'^^. 

2.1 DIS on a massive quark at 0{a^) 

The 0{a^) results for B*Qi Q2, including mass effects, have been obtained in ^ 
within a helicity basis for the hadronic/partonic structure functions. For completeness 
and in order to define our normalization we repeat these results here within the standard 
tensor basis implying the usual structure functions -Fi=i,2,3- The helicity basis seems to be 
advantageous since in the tensor basis partonic structure functions mix to give hadronic 
structure functions in the presence of masses . However, the mixing matrix is diagonal 
for the experimental relevant structure functions -Fj=i,2,3 and only mixes F4 with F5 
which are both suppressed by two powers of the lepton mass. We neglect target (nucleon) 
mass corrections which are important at larger values of Bjorken-x where heavy quark 
contributions are of minor importance. 

We consider DIS of the virtual Boson B* on the quark Qi with mass mi producing 
the quark Q2 with mass 1712- At order 0{a^) this proceeds through the parton model 
diagram in Fig. 1 (a). Finite mass corrections to the massless parton model expressions 
are taken into account by adopting the Ansatz given in Eq. (4) of 



W'^'^ is the usual hadronic tensor and o)^^ is its partonic analogue. Here as in the following 
a hat on partonic quantities refers to unsubtracted amplitudes, i.e. expressions which still 




(1) 
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contain mass singularities in the massless limit, pf and are the light-cone momentum 
components of the incident quark Qi and the nucleon, respectively. Generally the '+' 
light-cone component of a vector v is given by = (f ° + v^) / ^/2. 

Contracting the convolution in Eq. (|l]) with the projectors in Appendix A gives the 
individual hadronic structure functions -Fi=i,2,3- In leading order (LO) the latter are given 
by i 

F^'''\x,Q') = 2R^Q^{x,Q') (2) 

with 

E±± = Q^±ml±ml . (3) 

In Eq. (0) we use the shorthand A = A[mf,m2, —Q"^] , where the usual triangle function 
is defined by 

A[a, b, c] = v^a^ + b"^ + c"^ - 2{ab + bc + ca) . (4) 

The vector (V) and axial vector (A) couplings of the Q27m(^ ~ ^l5)Qi quark current 
enter via the following combinations: 

S± = VV'±AA' 

R± = {VA'±V'A)/2 (5) 

where V,A = V, A' in the case of pure B scattering and V,Ay^ V, A' in the case of 
B, B' interference (e.g. 7, interference in the standard model). The scaling variable x 
generalizes the usual Bjorken-x in the presence of parton masses and is given by [^: 
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The mass dependent structure functions in Eq. (|^) motivate the following definitions 

Fi 



^1 

^3 



2A 



2Q2 1 p 



> 



5+A 2x " ^ 
2iJ+ -^3 

such that JFj — jFj, j = 1, 2, 3, will be finite of 0{as) in the limit mi 2 — > 0. 



(7) 



2.2 DIS on a massive quark at 0{al] 



At 0{al) contributions from real gluon emission [Fig. 1 (b)] and virtual corrections [Fig. 
1 (c)] have to be added to the 0{a^) results of section 2.1. The vertex correction with 
general masses and couplings [Fig. 1 (c)] does to our knowledge not exist in the literature 
and is presented in some detail in Appendix B. The final result (virtual+real) can be cast 
into the following form: 



l2 ,.2\ 



27r 



with 



Hl{i\rm,m2) = CF 



{Si + V{) 5(1-0 + 



^ — i' s — ml 



8s 



(■I X 



(9) 



where s = (pi + g)^ and the 6*4, V^, Ni and j'^ are given in Appendix C. The factorization 
scale /i^ will be taken equal to the renormalization scale throughout. The '+' distribution 
in Eq. (^ is a remnant of the cancellation of the soft divergencies from the real and virtual 
contributions. It is defined as usual: 



(10) 



'0 ^0 

As indicated by the hat on Hf, the full massive convolution in Eq. (|^) still contains the 
mass singularity arising from quasi-collinear gluon emission from the initial state quark 



leg. The latter has to be removed by subtraction in such a way that in the asymptotic 
limit Q'^ ^ oo the well known massless MS expressions are recovered. The MS limit is 
mandatory since all modern parton distributions - and therefore all available heavy quark 
densities - are defined in this particular scheme (or in the DIS scheme [|18| , which can be 
straightforwardly derived from MS). The correct subtraction term can be obtained from 
the following limit 

X 



lim , 



+ C 



1-e 



1 -21n(l 



(11) 



where A = Q^/iQ^ + ml), x/X — x\mi=o and the Hf' can be found in Appendix A of 



p5| . Obviously the MS subtraction term for a 'heavy quark inside a heavy quark' is given 
not only by the splitting function Pgq^ = Cf[{^ + ^'^)/(l — ^')]+ times the collinear log 
ln(/i^/mf) but also comprises a constant term. Herein we agree with Eq. (3.15) in , 
where this was first pointed out in the framework of perturbative fragmentation functions 
for heavy quarks. We therefore define 



^^^^(a:,Q^/.^ 



such that 



1+e 



l2 



In 



- 1 -21n(l 



X 2 



lim 

mi— >0 



MS, 



x,Q^/i^) 



(12) 



(13) 



MS 



where the superscript Qi on JF. ^ ' refers to that part of the inclusive structure function 
J-'i which is initiated by the heavy quark Qi, i.e. which is obtained from a convolution 
with the heavy quark parton density. Note that the limit in Eq. (|TT|) guarantees that Eq. 
(13) is also fulfilled when mi = m2 — (e.g. NC leptoproduction of charm) since 



(14) 



^ We also agree with the quark initiated coefficient functions in |ll7| where quark; masses have been 
used as regulators 
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where Cf'^^ are the standard massless coefficient functions in the MS scheme, e.g. in 



2.2.1 Comparison to existing NC and CC results 

We have performed several cross checks of our results against well known calculations that 



exist in the literature |15, O, 18, 19, BQ]. The checks can be partly inferred from section 



2.2. Nevertheless we present here a systematic list for the reader's convenience and in 



order to hint at several errors which we uncovered in 20 



In the charged current case V = A = 1 our results in Eq. (B) reduce in the limit 



mi — i> to the corresponding expressions in [T^, or in [TH] if the scheme dependent term 



represented by Eq. (|T^) is taken into account. The latter agrees with Eq. (3.15) in ||T6| . 
For mi 2 ^ we reproduce the well known MS coefficient functions e.g. in [l^, [T^. The 
vertex correction in Appendix C is implicitly tested because it contributes to any of the 
final results. However, as an independent cross check the well known QED textbook result 
can be reproduced for mi = m2, A = 0. 



Initial state parton mass effects in NC DIS at 0{al) have been ffist considered in 



within the scenario of intrinsic nonperturbative cc pairs stemming from fluctuations of 
the nucleon Fock space wavefunction. Although we do not consider such a scenario here 
we note that our results could be easily transferred to corresponding applications |22 . 



The main difference would be an inclusion of kinematical target mass effects which are 
important at larger x where a possible nonperturbative charm component is expected 
pT| to reside. Apart from obvious typos, we uncovered some errors in such that our 
results differ from those in |^ by terms that vanish in the massless limit. Eq. (C.16) in 
should be multiplied by (1 + 4Az^). The typo propagates into the final result in Eq. 



(51). Furthermore the longitudinal cross section, i.e. Eq. (C.17) in pO[ seems to have been 
obtained as a residual of cr^^^ and the (wrong) a^^\ Constructing cr^'' from Eqs. (C.16), 
(C.17) via cr^'* = — 2(T^^ + crj^"* reproduces up to a constant of normalization the part 



7 



~ fi of our result for Fi in Eq. (|). Given the amount of successful independent tests 
of our results we regard the disagreement with ||20| as a clear evidence that the results in 
|20| should be updated by our calculation. 



2.3 Gluon fusion contributions at 0{al) 

The gluon fusion contributions to heavy quark structure functions {B*g —>■ Q1Q2) are 
known for a long time ||2^, ^ and have been reinterpreted in ||10| within the helicity 
basis for structure functions. Here we only briefly recall the corresponding formulae in 
the tensor basis for completeness. The GF component of DI structure functions is given 
by 

Frix,Q') = £ |^r^7(^',/^^)/2(^,Q^) (15) 

where ax = [1 + (mi +7/12)^/ Q'^]x and the fi can be found for general masses and couplings 
23|1 . The corresponding J-"^^ are obtained from the F[^^ by using the same normal- 



m 



ization factors as in Eq. (|^). Along the lines of [|I^ the GF contributions coexist with 
the QS contributions which are calculated from the heavy quark density, which is evolved 
via the massless RG equations in the MS scheme. As already pointed out in section 2.2 
the quasi-coUinear log of the fully massive GF term has to be subtracted since the cor- 
responding mass singularities are resummed to all orders in the massless RG evolution. 
The subtraction term for the GF contribution is given by ||10[ 

Q\ = E ^ £j /' f PSHO S (f . . (16) 

where Pi?(^') = 1/2 [C^ + (1 - C?]- Note that Eq. ([l|) as well as Eq. ([ID are defined 
relative to the jFj in Eq. (|^ and not with respect to the experimental structure functions 
Fi. The sum in Eq. ( |T6|) runs over the indices of the quarks Qk for which the quasi- 
collinear logs are resummed by massless evolution of a heavy quark density, i.e. k = 1, 
k = 2 OT k = 1,2. 
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2.4 ACOT structure functions at 0{al) 



As already mentioned in the introduction it is not quite clear how the perturbation series 
should be arranged for massive quarks, i.e. whether the counting is simply in powers 
of as as for light quarks or whether an intrinsic heavy quark density carries an extra 
power of as due to its prehistory as an extrinsic particle produced by pure GF. We are 
here interested in the QS^^^ component of heavy quark structure functions. Usually the 
latter is neglected in the ACOT formalism since it is assumed to be suppressed by one 
order of as with respect to the GF contribution as just explained above. It has however 
been demonstrated in [^5[ ^ within MS that this naive expectation is quantitatively not 
supported in the special case of semi- inclusive production of charm (dimuon events) in CC 
DIS. We therefore want to investigate the numerical relevance of the QS^^^ contribution 
to general heavy quark structure functions. In this article we present results for the fully 
inclusive case, relevant for inclusive analyses and fits to inclusive data. We postpone 
experimentally more relevant semi-inclusive (z-dependent) results to a future publication 
p7| . Our results at full 0{al) will be given by 

^(1) ^ ^Q5(o+i) _^ j^GF _ pSUB^ _ pSUBg ^-^^^ 

with F^^^°*'\ F^^, F-^^" and F^^^^' given in Eqs. (D, (|12D, and (0), respectively. 
Furthermore, we will also consider a perturbative expression for Fi which is constructed 
along the expectations of the original formulation of the ACOT scheme, i.e. QS^^^ is 
neglected and therefore pf^^'' need not be introduced 

p(0)+GF-SUBg ^ pQS(O) ^ pGF _ pSUB, ^^g^^ 



3 Results for NC and CC structure functions 

In this section we present results which clarify the numerical relevance of QS^^^ contribu- 
tions to inclusive heavy quark structure functions in the ACOT scheme. We will restrict 
ourselves to NC and CC production of charm since bottom contributions are insignificant 
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to present DI data. Our canonical parton distributions for the NC case will be CTEQ4M 
||28| (Figs. 4 and 5 below), which include 'massless heavy partons' Qk above the scale 
= ml. Figs. 2 and 3, however, have been obtained from the older GRV92 dis- 
tributions. The newer GRV94 parametrizations do not include a resummed charm 
density since they are constructed exclusively along FOPT. GRV94 is employed in the CC 
section. The radiative strange sea of GRV94 seems to be closest to presently available CC 



charm production data |2^. Furthermore, the low input scale of GRV94 allows for a wide 
range of variation of the factorization scale around the presently relevant experimental 
scales, which are lower for CC DIS than for NC DIS. Qualitatively all our results do not 
depend on the specific set of parton distributions chosen. 

3.1 NC structure functions 

For our qualitative analysis we are only considering photon exchange and we neglect the 
Z^. The relevant formulae are all given in section 2 with the following identifications: 



Q 



1,2 



mi,2 = 1.6 (1.5) GeV for CTEQ4 (GRV92) 

V = V',A = A' -,0 

3' 

and we use /i^ = if not otherwise noted. We consider contributions from charmed 
quarks and anti-quarks which are inseparably mixed by the GF contribution. This means 
that in Eq. (|16D the sum runs over k = 1,2 and the relevant expressions of section 2.1 
and 2.2 have to be doubled [since c{x,ij?) = c(x, /i^)]. 

First we investigate the importance of finite mass corrections to the limit in Eq. (|13D. 
In Fig. 2 the difference — F2 can be compared to its MS analogue which is 



(i),MS 4 as(/i2 



^ 9 271 



(c + c)(/x2)®C^ 



,MS / Q 



")2 



(x,g2) (19) 



where ® denotes the usual (massless) convolution. From Fig. 2 it is obvious that the 
relative difference between ACOT and MS depends crucially on x. It can be large and 
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only slowly convergent to the asymptotic MS limit as can be inferred from Fig. 3. Note 
that the solid curves in Figs. 2, 3 are extremely sensitive to the precise definition of the 
subtraction term in Eq. (|12D, e.g. changing x — > x - which also removes the collinear 
singularity in the high limit - can change the ACOT result by about a factor of 5 
around ~ 5 GeV^.0 This is an example of the ambiguities in defining a variable flavor 



number scheme which have been formulated in a systematic manner in |11 



The relative difference between the subtracted QS^^^ contribution calculated along 
ACOT and the corresponding MS contribution in Eq. (|l^) appears, however, phenomeno- 
logically irrelevant if one considers the significance of these contributions to the total 
charm structure function in Fig. 4. The complete 0[a^^ result (solid line) is shown over 
a wide range of together with its individual contributions from Eq. ([l^). It can be 
clearly seen that the full massive QS^^^ contribution is almost completely canceled by the 
subtraction term SUBq (Indeed the curves for QS^^^ and SUBq are hardly distinguish- 
able on the scale of Fig. 4) . The subtracted quark correction is numerically negligible and 
turns out to be indeed suppressed compared to the gluon initiated contribution, which is 
also shown in Fig. 4. Note, however, that the quark initiated corrections are not unimpor- 
tant because they are intrinsically small. Rather the large massive contribution QS^^^ is 
perfectly canceled by the subtraction term SU Bq provided that = is chosen. This 
is not necessarily the case for different choices of /x^ as we will now demonstrate. 

In Fig. 5 we show the dependence of the complete structure function and its compo- 
nents on the arbitrary factorization scale /i^. Apart from the canonical choice /i^ = 
(which was used for all preceding figures) also different scales have been proposed ^ 
like the maximum transverse momentum of the outgoing heavy quark which is approx- 
imately given by [p^"-^)^ ^ (1/a; — 1) Q^/4. For low values of x, where heavy quark 
structure functions are most important, the scale (p™"^)^ ^ Q^. The effect of choosing 
a yU^ which differs much from can be easiest understood for the massless coefficient 
functions C*?'^'^^ which contain an unresummed ln(Q^//i^). The latter is of course absent 
The subtracted gluon contribution GF changes by about a factor of 2 under the same replacement. 
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for /i^ = but becomes numerically increasingly important, the more /i^ deviates from 
Q^. This logarithmic contribution cannot be neglected since it is the unresummed part 
of the coUinear divergence which is necessary to define the scale dependence of the charm 
density. This expectation is confirmed by Fig. 5. For larger values of //^ the subtracted 
QS^^"* contribution is indeed still suppressed relative to the subtracted GF contribution. 
Nevertheless, its contribution to the total structure function becomes numerically signifi- 
cant and reaches the ~ 20 % level around (p^"^)^. Note that in this regime the involved 
formulae of section 2.2 may be safely approximated by the much simpler convolution in 
Eq. (^) because they are completely dominated by the universal collinear logarithm and 
the finite differences ACOT — MS from Figs. 2 and 3 become immaterial. In practice it is 
therefore always legitimate to approximate the ACOT results of section 2.2 by their MS 
analogues because both are either numerically insubstantial or logarithmically dominated. 
Finally we confirm that the scale dependence of the full 0{a\) structure function f!^^ 
in Eq. (|l^ is larger than the scale dependence of p^^)'^^^~^^^3 \yi Eq. (|1^) which was 
already pointed out in [^. Nevertheless, the subtracted QS'^^'> contribution should be 



respected for theoretical reasons whenever asln(Q^//U^) ^ 1. 
3.2 CC structure functions 

Charm production in CC DIS is induced by an s — >■ c transition at the VT-Boson vertex. 
The strange quark is not really a heavy quark in the sense of the ACOT formalism, i.e., 
the production of strange quarks cannot be calculated reliably at any scale using FOPT 
because the strange quark mass is too small. It is nevertheless reasonable to take into 
account possible finite effects into perturbative calculations using ACOT since the 
subtraction terms remove all long distance physics from the coefficient functions. Indeed 
the ACOT formalism has been used for an experimental analysis of CC charm production 
in order to extract the strange sea density of the nucleon |^. Along the assumptions of 
ACOT QS'^^^ contributions have not been taken into account. This procedure is obvi- 
ously questionable and has been shown not to be justified within the MS scheme |2^ . 



12 



With our results in section 2.2 we can investigate the importance of quark initiated 0{al) 
corrections within the ACOT scheme for inclusive CC DIS. As already mentioned above, 
results for the experimentally more important case of semi-inclusive (2;-dependent) DIS 
will be presented in a future publication In the following we only introduce sub- 



traction terms for collinear divergencies correlated with the strange mass and treat all 
logarithms of the charm mass along FOPT. We do so for two reasons, one theoretical 
and one experimental: First, at present experimental scales of CC charm production 
\n{Q'^/ml) terms can be safely treated along FOPT and no introduction of an a priori 
unknown charm density is necessary. Second, the introduction of a subtraction term for 
the mass singularity of the charm quark would simultaneously require the inclusion of 
the c — s> s QS'-transition at the W-vertex with no spectator-like c-quark as in GF. This 
contribution must, however, be absent when experiments tag on charm in the final state. 
CC DIS on massive charm quarks without final-state charm tagging has been studied in 



m 



The numerics of this section can be obtained by the formulae of section 2 with the 
following identifications: 

Qi ^ s , Q2 ^ c 
m2^mc = 1.5 GeV (GRV94) 
V = V',A = A' 1,1 

and the strange mass mi = will be varied in order to show its effect on the structure 
function F2. 

In Fig. 6 we show the structure function F| and its individual contributions for two 
experimental values of x and under variation of the factorization scale fi^. Like 
in the NC case we show the complete 0{al) result as well as ^2^^'^'^^ suBg ^^lexe QS^^'' 
has been neglected. The thick curves in Fig. 6 (a) have been obtained with a regularizing 
strange mass of 10 MeV. They are numerically indistinguishable from the = MS 
results along the lines of [|1^. For the thin curves a larger strange mass of 500 MeV 
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has been assumed as an upper limit. Finite mass effects can therefore be inferred from 
the difference between the thin and the thick curves. Obviously they are very small for 
all contributions and can be safely neglected. For the higher value of Fig. 6 (b) they 
would be completely invisible, so we only show the =10 MeV results (=MS). Since the 
finite mass corrections within the ACOT scheme turn out to be negligible as compared 
to massless MS it is not surprising that we confirm the findings of p5| , |26| concerning 
the importance of quark initiated corrections. They are - in the case of CC production 
of charm - not suppressed with respect to gluon initiated corrections for all reasonable 
values of the factorization scale. Only for small choices of /i^ ~ + '^^'^ quark 
initiated correction be neglected. In this region of //^ also gluon initiated corrections are 



moderate and Born approximation holds within ~ 10% |T^. For reasons explained in 
section 3.1 the absolute value of both corrections - gluon and quark initiated - become 
very significant when large factorization scales like p^"^ are chosen. This can be inferred 
by looking at the region indicated by the arrow in Fig. 6 which marks the scale // = 2p^°'^ 
which was used in 0. Analyses which use ACOT with a high factorization scale and 
neglect quark initiated corrections therefore undershoot the complete 0{al) result by the 
difference between the solid and the dot-dashed curve, which can be easily as large as ~ 20 
%. For reasons explained in the introduction to this section we have used the radiative 
strange sea of GRV94. When larger strange seas like CTEQ4 are used the inclusion of 
the quark initiated contributions is even more important. 



4 Conclusions 



In this article we have calculated and analysed DIS on massive quarks at 0{al) within the 
ACOT scheme for heavy quarks. For NC DIS this contribution differs significantly from its 
massless MS analogue for /x^ = Q^. Both give, however, a very small contribution to the 
total charm structure function such that the large relative difference is phenomenologically 
immaterial. At higher values of the factorization scale /i^ ~ (p™'^^)^ the contributions 
become significant and their relative difference vanishes. The QS^^^ contribution of section 
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2.2 can therefore be safely approximated by its much simpler MS analogue at any scale. 
For CC DIS quark initiated corrections should always be taken into account on the same 
level as gluon initiated corrections. Due to the smallness of the strange quark mass ACOT 
gives results which are almost identical to MS. 
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Appendix A: Real Gluon Emission 



We define a partonic tensor 



color spin 

which can be decomposed into its different tensor components as usual 



-Cor 9'" + ^7 P'^Pi + ejp^.q^ + q'q" + Co. 



(Al) 



(A2) 



Co^y can be easily calculated from the general Feynman rules for invariant matrix elements 
which are customarily expressed as functions of the mandelstam variables s = (pi + g)^ 
and i = {pi — kY to which we will refer in the following. Projection onto the individual 



00. 



Q 

i=l,2,3 



in Eq. 



is performed for nonzero masses and in = 4 + 2e dimensions with 



the following operators 
-1 



)2 H V 



-'2 



2(1+^; 

- g^'^'Q^ + 4 q^q 



, 2(l+e)(pi.g)2-m2g2 



^3 



+ A{3 + 2e)Q 
-2i 



;i + £)A'K,s,-g 



«i ^ xk Pi q 



(A3) 



A2K,s,-g2 

such that Pi ■ Co = CJ^ . The normalization in Eqs. ( [A2|) is such that real gluon 

emission contributes Ff- to the hadronic structure functions via 

1 /"^ di 



R 



Fi 



2x di A2[m2,s,-Q2 

7 



IGtt 
1 

8^ 



Qi(0 / rfPS ^2^2^ 



2g2 

AK,S,-Q'] Qi(0 / ^^PS^3 



where 35 



c/PS = — 



1 s — m\ 



s V{l + e) 



{s — m. 
47rs 



2\2 



(A4) 



(A5) 
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is the partonic phase space. In Eq. ( |A5| ) y is related to the partonic centre of mass 
scattering angle 6* and the partonic mandelstam variable i via 



y = \ (i + cosr 



rrin 



(A6) 



We have chosen dimensional regularization for the soft gluon poles stemming from s 



which arise from propagators in the Ui times phase space factors in dPS. In Eq. (|A4|) we 
use 



18 



(s — m 



2\2e-l 



1 



2e-l 



which separates hard gluon emission (~ f^) from soft gluon (Si) contributions in Eq. (^. 
Note that in Eqs. (|]), (|^) the integration variable ^, which is implicitly defined in Eq. (|l]), 
has been changed to C,' = x/^ fo^' easier handling of the distributions. For the relation 
between s and ^' see Eq. ( p3|) below. 



Since all quark masses are kept nonzero, no poles in y (collinear singularities) are 
contained in the integration volume. The which occur in Eq. (^) and which are given 
below in Eq. ( |C4]) are therefore straightforward integrals of the 

f?={9lC,)-' f dyu^ (A8) 
Jo 

and the Si in Eqs. (1), ( PI]) pick up the pole in Eq. ( |A^ ) 

1 



Si 



dy [y{l-y)Y Cj)^{s-mlf 



(A9) 



where the proportionality is given by kinematical and phase space factors which must be 
kept up to 0{e). 



The normalization of our hadronic structure functions in Eq. (|A4|) can be clearly 

inferred from the corresponding LO results in Eq. (Q). Nevertheless, for definiteness we 

also give the hadronic differential cross section to which it corresponds 

^2-''^ l_ {Gf'^^'f {Gf^'f 2M^E, 
ni 



d^a' 



dxdy 



27r 



Si,+ (1 - y)F2 + Si,+ y'xF^ ± 2y(l - ^)xFs 



(AlO) 
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where {Gfff = [<V(Q^ + ^s) 9?'^ I{Q'^ + ^I')]^^^ is the effective squared gauge 
couphng - including the gauge boson propagator - of the 7^(1^ — ^475) lepton and quark 
current, respectively, and ni counts the spin degrees of freedom of the lepton, e.g. = 1, 2 
for I = v,e~ . The leptonic couplings 5; +, i?/ + are defined analogous to the quark couplings 
in Eq. (|^). As noted below Eq. (^), B = B' for non-interference (pure B scattering). 



Appendix B: Vertex Correction 

We have calculated the vertex correction in n = 4 + 25 dimensions at 0{a\) for gen- 
eral masses and couplings using the Feynman gauge. The unrenormalized vertex Aq 
[Fig. l(c.l)] has the structure 



AM 
^^0 



a 



Co,- i^L, + C+ 7'^^5 



+ m2 L5 + mi i?5 + Cg,_ msg^ L5 + Cg,+ rmg^ i^sj (Bl) 



with L5 

Co,- ^ 



(V^ — A 75), i?5 = (V^ + A 75). The coefficients read 



+ 



X 



Ci,- 



-(-1 - S^+Ji) 



2g2 



2 \m{j 2 



-In- 
2 

Li2 



A - 



2g2 

A-S+_ 



+ -ln^ 



2g2 

A - S_ 



m? / m 
- m 

mi 



h 



2A 



Li2 



2Q^ 
A - 
2A~ 



-In' 



1 \ _|_ ^++ 



A + S, 



(B2) 



2 

+ Li2 



2g2 

A-hS+_ 
2A 



-In^ 



A + S_, 



+ Lis 



2g2 

A + S_, 
2A 



2mim2/i 
-1 



Cn4- — — — r 



-1 

g^ 

1 

g^ 
1 



S,_/i + In f 4 



S_+Ji - In 



/ 2 
/ mi 



(A^ - 2m2g2) - 2g2 + In ^ 

(-A^ + 2m2g2 - s_+g2) ji + 2g2 + (s_ 



(B3) 
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with 



h = JilnA- — 
A 



X 



2 
Li2 



In- 



A-E, 



A-E+_ 
2A 



A-E_ 



Lio 



2Q2 
A-E_+ 
2A 



2g2 



2 



L12 



A + E, 



2g2 



V 2A 



2 



T- /A + E_ 
L12 



A + E_ 



V 2A 



(B4) 



The renormahzed vertex [Fig. 1 (c.l)-(c.3)] is obtained by wave function renormahza- 
tion: 

K^^^YL,{Z,-l)-rK + 0{al) (B5) 



where Zi = a/^2(Pi)^2(P2)- The fermion wave function renormahzation constants are 
defined on mass sheU: 



such that 



2 \ £ 



3 3 /Q2 
in 



£ 2 



rrio 



(B6) 



(B7) 



The final result for the renormahzed vertex reads 

12 \ e 



+ Ci,_ m2 V^i L5 + Ci,+ mi p1 R5 + Cg,_ m2 q" L5 + Cg,+ mi q" i^gj (B8) 



with C+, Ci,±, Cg,± as given above and 



-(2-E++7i) + 



2g^ + ^-l^ + ^nl 



+ 



X 



m2 — mj 



2Q2 



m^ 



mt 



h 



4 + 



(B9) 



-In^ 



A-E, 



— Li9 



2g2 

A-E+. 
2A 



+ -ln2 
2 



A-E_+ 




A + E+_ 




A + E_+ 


2g2 


2 


2g2 


2 


2g2 



- Lio 



A-E_ 
2A~ 



+ Lic 



/A + E, 
V 2A 



+ Li. 



/A + E_ 



V 2A 
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Appendix C: Real and Virtual Contributions to Struc- 
ture Functions 



The soft real contributions Si to the coefficient functions in Eq. 

2A 



Si 



A Ji + Lis 



A-E, 



are given by 
2A 



Li-; 



A + S, 



S 



2,3 



^1 



(CI) 



with Ii given in Appendix B and where x is given in Eq. (^. The virtual contributions 
are derived from the renormalized vertex in Eq. ( pSj) by using the projectors in Eq. (|A3|) : 



V3 



- 2S+mim2 

— 2S-mim2 

Cr,- + ^ {ml + ml Ci,_) + |- 



R 



(C2) 



+ 



where the C's are given in Appendix B. Note that the soft poles of Si, Vi cancel in 

the sum Si + Vi in Eq. (|) as must be. 



The massive matrix elements f^{^') are most conveniently given as functions of the 
mandelstam variable = {pi + q)^ — ra\, i.e. f^{C) = f^l^iiO] with 

.2 1-^', 



s — rrir, 



2i 



^[(A-S^_)^' + A + S, 



(C3) 



From the real graphs of Fig. 1 (b) one obtains 



A'2 

+ 2E+. -E_+ + 



A^(5+S++ - 2mim2^_)/e + 2mim2^_ ( ^[A'^ + 4m^S^ 



+ si , si + m\ 



+ 



A'si 



+ ^4 



A'' + 2S+_S++ + [ml + 
1 



4(si +m2) 



- - m2s++) + si[-7S++S_+ + 18A2 - 4m2(7Q2 - Ami + 
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with 



and 



+ 
+ 



+ 



2mi 



Si + 1712 

2A' 



Sl 



A^"(gi + 

2(si + ml 
2 (m^ + ml) sl - 9m^E^_ + A^ (2S 



;) + 2si (2A2 + (m^ - 5ml) 



A'' - + Sl)) (si + 2A2 



2(si + ml) 



^^^+"^^)[Z^A^(A^ + 2S,_S, 



Sl 
.2 



(A2 + 2(2m2 + si)S+_) 



2si(A^ - 6m2Q2^ 



A' 



- (A'^ - 18m^Q")S++ - 2A" + 2S+_)] L^, 



^1 - 2A^R+I^, + 2mim2R- (^1 



+ 



+ i?+ ( S_+ - 3S+_ - ^ (A^ + 2m^E+_) 



Sl A'si 
Sl - E_+)(si + 



2(si + ml 



Sl + rrij 
A'si 



[-s? + 4(m2s_+- A^) -3siS+_] L^/ 



Lt/ = In 



Sl - A' 



E++ + Sl + A' 
Sl + 2ml Sl + ml 



si A'si 
A is given below Eq. (H) and A' = A[mi, s, —Q"^]. 



2 S++ -^e 



(C4) 



Finally, the normalization factors in Eq. 

S+T,++ — 2mim2S_ 



are 



iVi 



2A 



25+A 
(A')2 ^ 



AT, 



2i?^ 



(C5) 



21 



References 



[1] C. Adloff et ai, HI coUab., Z. Phys. C72, 593 (1996). 

[2] J. Breitweg et ai, ZEUS coUab., Phys. Lett. B407, 402 (1997). 

[3] W. G. Seligman et ai, Phys. Rev. Lett. 79, 1213 (1997); 

W. G. Sehgman, Ph. D. thesis, Columbia University, Nevis-292 (1997). 

[4] A. Vogt, proceedings of the International Workshop on Deep Inelastic Scattering and 
Related Phenomena (DIS 96), Rome, Italy, April 1996; Edited by G. D'Agostini and 
A. Nigro, World Scientific, 1997. 

[5] H. Abramowicz et ai, CDHSW coUab., C. Phys. C15, 19 (1982). 

[6] S. A. Rabinowitz et ai, CCFR coUab., Phys. Rev. Lett. 70, 134 (1993). 

[7] A. O. Bazarko et ai, CCFR coUab., Z. Phys. C65, 189 (1995); 

A. O. Bazarko, Ph. D. thesis, Columbia University, Nevis-285 (1994). 

[8] M. Gliick, E. Reya and M. Stratmann, Nucl. Phys. B422, 37 (1994). 

[9] M. A. G. Aivazis, F. I. Olness and W.- K. Tung, Phys. Rev. D50, 3085 (1994). 

[10] M. A. G. Aivazis, J. C. Collins, F. I. Olness and W.- K. Tung, Phys. Rev. D50, 3102 
(1994); 

J. C. Collins PSU-TH/198, hep-ph 9806259. 

[11] R. S. Thorne and R. G. Roberts, RAL-TR-97-049, hep-ph 9709442; 
R. S. Thorne and R. G. Roberts, RAL-TR-97-061, hep-ph 9711223. 

[12] A. D. Martin, R. G. Roberts, M. G. Ryskin and W. J. Stirling, Eur. Phys. J. C2, 
287 (1998). 

[13] M. Buza, Y. Matiounine, J. Smith, W. L. van Neerven, Eur. Phys. J. CI, 301 (1998). 
[14] J. C. Collins and W.- K. Tung, Nucl. Phys. B278, 934 (1986). 

22 



[15] M. Gliick, S. Kretzer and E. Reya, Phys. Lett. B380, 171 (1996); Erratum B405, 
391 (1996). 



[16: 
[17: 
[18 

[19 

[2o: 

[21 
[22 
[23: 
[24 

[25: 



[27 

[28: 
[29 

[3o: 



B. Mele and P. Nason, Nucl. Phys. B361, 626 (1991). 

J. J. van der Bij and G. J. Oldenborgh, Z. Phys. C51, 477 (1991). 

G. Altarelh, R. K. Elhs and G. Martinelh, Nucl. Phys. B157, 461 (1979). 

W. Furmanski and R. Petronzio, Z. Phys. Cll, 293 (1982). 

E. Hoffmann and R. Moore, Z. Phys. C20, 71 (1983). 

S. J. Brodsky P. Hoyer, C. Peterson and N. Sakai, Phys. Lett. B93, 451 (1980); 
S. J. Brodsky, C. Peterson and N. Sakai, Phys. Rev. D23, 2745 (1981). 

B. W. Harris, J. Smith, R. Vogt, Nucl. Phys. B461, 181 (1996); 

J. F. Gunion and R. Vogt, UCD-97-14, LBNL-40399, hep-ph 9706252. 

M. Gliick, R. M. Godbole and E. Reya, Z. Phys. C38, 441 (1988); Erratum C39, 
590 (1988). 

G. Schuler, Nucl. Phys. B299, 21 (1988); 

U. Baur and J. J. van der Bij, Nucl. Phys. B304, 451 (1988). 

M. Gliick, S. Kretzer and E. Reya, Phys. Lett. B398, 381 (1997); Erratum B405, 
392 (1997). 

S. Kretzer and I. Schienbein, Phys. Rev. D56, 1804 (1997). 
S. Kretzer and I. Schienbein, (in preparation). 

H. L. Lai et al, CTEQ coUab., Phys. Rev. D55, 1280 (1997). 
M. Gliick, E. Reya and A. Vogt, Z. Phys. C53, 127 (1992). 
M. Gliick, E. Reya and A. Vogt, Z. Phys. C67, 433 (1995). 



23 



[31] M. A. G. Aivazis, F. I. Olness and W.- K. Tung, Phys. Rev. Lett. 65, 2339 (1990). 

[32] J. C. Collins, Proceedings of the XXVth Rencontre de Moriond: High Energy 
Hadronic Interactions, ed. J. Tran Than Van, Gif-sur-Yvette: Editions Prontieres, 
1990, page 123. 

[33] C. R. Schmidt, proceedings of the 5th International Workshop on Deep Inelastic 
Scattering and QCD (DIS 97), Chicago, IL, USA, April 1997. 

[34] U. D'Alesio, Ph. D. thesis {in Italian), Turin University (1996). 

[35] T. Gottschalk, Phys. Rev. D23, 56 (1981). 



24 



Figure Captions 



Fig. 1 Feynman diagrams for the Q^(o) [Fig. 1 (a)] and QS^^^ [Fig. 1 (b), (c)] contribu- 
tions to ACOT structure functions in Eqs. @ and (H), respectively. 

Fig. 2 x-dependence of the subtracted QS^^^ contribution to the NC charm structure 
function (sohd hne). Q"^ = fi^ = 10 GeV^ is fixed. For comparison the MS 
analogue in Eq. ( |T9| ) is shown (dashed line). The GRV92 parton distributions have 
been used. 

Fig. 3 The same as Fig. 2 but varying Q^(= yU^) for fixed x. 

Fig. 4 The complete 0{al) neutral current structure function and all individual con- 
tributions over a wide range of Q^, calculated from the CTEQ4M distributions. 
Details of the distinct contributions are given in the text. 

Fig. 5 /i^ dependence of the complete 0{al) NC structure function (CTEQ4M) in Eq. 
([TTp (solid line) and of the structure function in Eq. (|T^) (dot-dashed line) where the 
subtracted QS^^^ contribution is neglected. Also shown are the different subtracted 
0{al) contributions GF and QS^^\ 

Fig. 6 (a), (b) The charm production contribution to the charged current structure 
function F2 for a wide range of the factorization scale fi^ using GRV94. The curves 
are as for the neutral current case in Fig. 5. In Fig. 6 (a) the thicker curves have 
been obtained with a (purely regularizing) strange mass of 10 MeV which according 
to Eq. (0) (and to the analogous limit for the subtracted GF term |jTO|) numeri- 
cally reproduces MS. For the thinner curves a strange mass of 500 MeV has been 
assumed. In Fig. 6 (b) all curves correspond to = 10 MeV (= MS). 
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